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ABSTRACT 



VirB activates transcription of virulence genes 
in Shigella flexneri by alleviating heat-stable 
nucleoid-structuring protein-mediated promoter 
repression. VirB is unrelated to the conventional 
transcriptional regulators, but homologous to the 
plasmid partitioning proteins. We determined the 
crystal structures of VirB HTH domain bound by 
the c/s-acting site containing the inverted repeat, 
revealing that the VirB-DNA complex is related to 
ParB-ParS-like complexes, presenting an example 
that a ParB-like protein acts exclusively in transcrip- 
tional regulation. The HTH domain of VirB docks 
DNA major groove and provides multiple contacts 
to backbone and bases, in which the only specific 
base readout is mediated by R167. VirB only recog- 
nizes one half site of the inverted repeats containing 
the most matches to the consensus for VirB binding. 
The binding of VirB induces DNA conformational 
changes and introduces a bend at an invariant 
A-tract segment in the c/s-acting site, suggesting a 
role of DNA remodeling. VirB exhibits positive 
cooperativity in DNA binding that is contributed by 
the C-terminal domain facilitating VirB oligomeriza- 
tion. The isolated HTH domain only confers partial 
DNA specificity. Additional determinants for se- 
quence specificity may reside in N- or C-terminal 
domains. Collectively, our findings support and 
extend a previously proposed model for relieving 
heat-stable nucleoid-structuring protein-mediated 
repression by VirB. 



INTRODUCTION 

The pathogenic bacterium evolves from its harmless 
ancestor through the acquisition of the virulence genes 
located on plasmids, bacteriophages or pathogenic 
islands (1). The expression of the virulence genes enables 
the microbe to survive in new host environment. The evo- 
lution of Shigella from the innocuous Escherichia coli rep- 
resents a good example for the understanding of how the 
acquired virulence genes contribute to pathogenic lifestyle 
of the bacterium. Shigella is the causative agent for bacil- 
lary dysentery, posing a leading public health concern 
worldwide (2,3). Shigella can be subdivided into four 
species, Shigella flexneri, Shigella dysenteriae, Shigella 
sonnei and Shigella boydii, in which most of the pathogen- 
icity studies were carried out on the prevalent S. flexneri. 
The virulence of S. flexneri depends on the virulence genes 
clustered within the so called 'entry region', a 31-kb 
segment on the virulence plasmid (4,5). These genes 
encode a variety of proteins, including invasins mediating 
the invasion of intestinal epithelium, structural proteins 
forming the type III secretion system and the effector 
proteins (6). The vast quantity of the virulence genes 
products presents a great metabolic burden; therefore, 
their expressions are strictly controlled involving both 
the plasmid encoded regulatory system and the chromo- 
some encoded global gene regulation network. When 
the growth conditions are not ideal for invasion, the tran- 
scription of the virulence genes is repressed by a 
chromosomally encoded heat-stable nucleoid-structuring 
protein (H-NS) (7). H-NS binds the AT-rich DNA 
segment near the promoters, from which H-NS spreads 
along the DNA, forming H-NS-DNA filaments or 
DNA-H-NS-DNA bridges (8,9). The resulting nucleopro- 
tein structures function as the obstacles to impede the 
movements of RNA polymerase, thereby silent the 
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transcription (10,1 1). In response to proper environmental 
signals, i.e. pH value, temperature and osmolarity similar 
to the site of invasion (12), the H-NS-mediated transcrip- 
tional repression can be alleviated. A transcriptional 
cascade is then initiated with the activation of virF gene 
expressing an AraC-like protein VirF that in turn activates 
the transcription of the virB gene (13-15). The product 
VirB relieves the H-NS-mediated transcriptional repres- 
sion, leading to the activation of many virulence genes 
(16). 

VirB belongs to an unusual type of transcriptional ac- 
tivator because the protein does not recruit RNA poly- 
merase, indicating an indirect transcriptional activation 
(17,18). Indeed, VirB shares little sequence homology 
with the conventional transcriptional regulators, but it 
is homologous to the plasmid partitioning proteins 
ParB, SopB and KorB (19-21). VirB-mediated transcrip- 
tional activation initiates with the recognition of the ex- 
acting sites upstream promoter. It was speculated that on 
the binding to the c/s-acting site, VirB could destabilize 
and collapse the repressive H-NS-DNA complexes, 
releasing of the trapped RNA polymerase so that the 
transcription reactivates (17,22,23). However, the mech- 
anism underlying the VirB-mediated promoter derepres- 
sion remained elusive. Previous studies identified several 
c/s-acting sites for VirB binding. A consensus [5 / -(A/ 
G)(A/T)G(G)AAAT-3 / ] was derived from the analyses 
of InvE (identical to VirB) binding sites upstream icsB, 
spal5 and virA promoters in S. sonnei (24) and the VirB 
binding sites upstream the icsB, icsP and virB promoters 
in S. flexneri (17,23,25). The consensus shares the core 
elements with the A-boxes in ParS-like site targeted by 
the ParB-like proteins (21,26). Especially, the c/s-acting 
sites upstream icsB and icsP promoters contain the 
inverted repeats similar to the A2-A3 boxes in ParS site 
(21). 

The investigation of VirB-DNA interaction is the key 
to elucidate of the function of VirB. Similar to many 
sequence specific DNA-binding proteins, VirB binds 
DNA nonspecifically in vitro (22,27), whereas the recogni- 
tion of the specific DNA sequence by VirB is essential for 
transcriptional activation in vivo (24). Mutations at the 
VirB binding site could lead to the complete loss of 
VirB-mediated promoter activation (17,25). However, 
the structural basis for the DNA sequence specificity is 
unavailable. 



MATERIALS AND METHODS 

Construction of plasmids 

To construct plasmids encoding VirB variants and 
mutants, the XhoI-EcoRI fragment was generated by 
PCR and inserted into pBAD-mychis (Invitrogen). Point 
mutations were introduced using site-directed mutagenesis 
(QuickChange™). The icsP promoter (a 1256-bp 
fragment upstream of the icsP transcription start site) 
was amplified from the virulence plasmid pCP301 and 
combined to a promoterless lacZ gene and inserted into 
the lower-copy-number plasmid pACYC184 with BamHI 



and Xbal site. The sequences were verified by DNA 
sequencing. 

Gene mutation by lambda red recombination 

The virB gene knockout mutant strain was made by hom- 
ologous recombination using the A-red recombination 
system described previously (28). Briefly, the kanamycin 
resistance gene from pKD13 plasmid was amplified by 
PCR using primer pairs whose sequences are homologous 
to the 45 bp regions upstream and downstream of the virB 
gene. The PCR products were purified and transformed 
into electro-competent S. flexneri 2a strain 301 containing 
the red helper plasmid pKD46. The structure of virBr.kan 
lesion was confirmed by PCR and DNA sequencing. 

P-galactosidase assay 

Transcription of the LacZ gene under the control of icsP 
was determined by measuring P-galactosidase activity in 
bacteria S. flexneri 2a virB knockout strain containing 
plasmids pACYC184 and pBAD-mychis expressing VirB 
variants and mutants as described previously (22). Assays 
were performed at least in triplicate, and the data are ex- 
pressed as the mean of two independent measurements. 

Protein expression, purification and crystallization 

Plasmids encoding VirB variants VirB full length (FL) 
(1-309 aa), VirB delta N (129-309 aa), VirB delta C 
(1-250 aa) and VirB core (129-250 aa) or VirB mutants 
was transformed to Rosetta™ (DE3) competent cells 
(Novagen), respectively. The expression of the recombin- 
ant protein was induced by adding 0.5 mM IPTG when 
the culture reached OD 600 ~1.0. The selenomethionine 
derivatives were obtained by inducing protein expression 
in E. coli strain B834 (DE3) grown in LeM ASTER 
medium containing L-selenomethionine. The harvested 
bacteria cells were resuspended in lysis buffer containing 
50 mM Tris-HCl (pH 7.0), 2M NaCl and 30 mM imid- 
azole and disrupted by ultrasonication. The protein puri- 
fication is composed of three steps. The crude lysate were 
first loaded to Ni-NTA resin (invitrogen) and eluted with 
200 mM imidazole. The eluate was further purified using 
Hitrap Q HP 5 ml of column (GE healthcase) with the 
linear gradient of 150-1000 mM NaCl. The polishing 
step of the purification was the size exclusion chromatog- 
raphy using Superdex 75 HR 10/30 column (GE 
Healthcare) equilibrated with 20 mM Tris-HCl (pH 7.0), 
150mM NaCl and 2mM DTT. To crystallize VirB core/ 
DNA complexes, multiple oligonucleotides (23-33 bp) 
containing the c/s-acting sites identified upstream icsP or 
icsB promoters were synthesized, annealed with their com- 
plementary stands and then mixed with VirB core 
(12mg/ml) at molar ratio DNA:protein = 2:1, respect- 
ively. The mixtures were incubated on ice for 2h to 
allow the assembly of the nucleoprotein complexes 
before the crystallization trials. Only icsB-wt and icsP-wt 
(Supplementary Table SI) co-crystallized with VirB core. 
The crystal containing Set-Met labeled VirB core and 
5-bromouracil (5-BRU) incorporated icsB-wt was 
obtained by mixing 0.8 ul of sample with 1 jlxI of reservoir 
buffer containing 0.1 M MgCl 2 , 0.1 M Sodium Acetate 
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pH = 4.5, 22% PEG400, 3% dimethyl sulfoxide (v/v) in a 
vapor diffusion system at 22° C. The crystals containing 
native VirB core/icsB-wt and VirB core/icsP-wt were 
obtained using the similar approaches, except the reservoir 
buffer was optimized to 0.1 M MgCl 2 , 0.1 M Sodium 
Acetate pH 4.3, 30% PEG400, 8% 1,3-Butanediol (v/v) 
and 0.1 M MgCl 2 0.1 M Sodium acetate pH 4.5, 16% 
PEG400, respectively. The cryocooling of the crystals 
were carried out by the soaking in reservoir buffer con- 
taining 10% ethylene glycol 30-60 s before the flash 
freezing in liquid nitrogen. 

Structure determination 

Highly redundant single anomalous dispersion data to 
2.42 A on crystals containing Se-Met labeled VirB core/ 
5-BRU incorporated icsB-wt and data to 2.78 A on the 
crystals containing native VirB core/icsP-wt were collected 
at beamline X06DA SLS PSI (Villigen, Switzerland). Data 
to 2.43 A on crystals containing native VirB core/icsB-wt 
were collected at beamline BL17U at SSRF (Shanghai, 
China). The data sets were integrated and scaled using 
XDS package. The space group identified for icsB-wt con- 
taining crystals was ¥2{l{l\ for icsP-wt containing crystals 
was C2. One VirB core and one dsDNA (icsB-wt or icsP- 
wt) are present in the asymmetric unit. Programs SHLEX 
C/D/E were used to locate the heavy atoms (one selenium 
on VirB core and two bromines on icsB-wt) and to calcu- 
late the initial phases, leading to an interpretable electron 
density map. The manual model building was carried out 
by using the program Coot. The structures of the native 
VirB core/icsB-wt and VirB core/icsP-wt were solved by 
molecular replacement (Phaser, CCP4 package) using 
the structures of Se-Met labeled VirB core and 5-BRU 
incorporated icsB-wt as the searching models. 
Refinement of all structures were performed by using 
PHENIX, which included bulk solvent correction, individ- 
ual B factor refinement, simulated annealing and anomal- 
ous f and F refinement. Data collection and model 
statistics are summarized in Table 1. 

Fluorescence and anisotropy measurement 

Fluorescence anisotropy experiments were performed with 
a Fluoromax-4 fluorimeter (HORIBA Jobin Yvon), 
equipped with a Glan-Thompson prism polarizer. 
Briefly, 1.4ml of buffer [50 mM Tris-HCl (pH 7.5), 
150mM NaCl] containing 4nM fluorescently labeled 
DNA FAM-icsB-wt-50 bp (SangonBiotech Shanghai 
Co., Ltd, Supplementary Table SI) was pre-equilibrated 
in a quartz cuvette at 15°C;. Protein samples were added 
stepwise. After 3min equilibration, the anisotropy data 
were collected using an excitation wavelength of 493 nm 
and monitoring the emission at 518 nm. The band pass 
was 3nm for excitation and 3nm for emission. A 
maximum number of 10 trials were performed until 
<2% deviation of the signal was reached. 

Electrophoretic mobility shift assays 

To prepare the probe for electrophoretic mobility shift 
assay (EMSA), the synthetic DNA icsB-wt was 5'-labeled 
using [ 32 P]y- ATP (5000 Ci/mmol, Furui Biotech, Beijing, 



China) and purified using G-25 MicroSpin columns (GE 
healthcare). VirB variants (2uM), labeled icsB-wt (2nM), 
unlabeled competitor DNA (icsB-wt or icsB-mut supplied 
in 2x, 4x, 8x, 16x, 32x and 64x excess of VirB variant) 
were mixed. The mixtures were incubated at room tem- 
perature in the binding buffer containing 20 mM Tris-HCl 
(pH 8.0), 50 mM NaCl, 5mM MgCl 2 , 5% glycerol (w/v), 
lOng/ml poly dl-dC and 0.1 mg/ml bovine serum albumin 
for 30min. The resulting samples were resolved on 6% 
native PAGE in 0.5 x TBE buffer at 20V/cm for M h. 
The gels were visualized with Typhoon Trio Variable 
Mode Imager (GE healthcare). 

Molecular weight determination by size-exclusion 
chromatography 

Superdex 75 10/300GL column (GE healthcare) was 
equilibrated with the buffer containing 50 mM Tris-HCl 
(pH 8.0) and lOOmM NaCl, and calibrated using molecu- 
lar weight standards, y-globulin (158kDa), ovalbumin 
(45kDa), myoglobin (17kDa) and vitamin B12 
(1.35 kDa). The purified VirB variants (~1 mg/ml) were 
loaded to the column running at the flow rate of 
0.15ml/min. 

RESULTS 

The mode of DNA binding by VirB 

To investigate the molecular determinants underlying the 
DNA binding, we analyzed the DNA binding properties 
of VirB and a set of truncations (Figure 1A). We tested the 
full-length protein, denoted 'VirB FL', the truncation 
lacking C-terminal (1-250 aa), denoted 'VirB delta C, 
the truncation lacking N-terminal (129-309 aa), denoted 
'VirB delta N' and truncation lacking both C- and N- 
terminal (129-250 aa), denoted 'VirB Core'. All VirB 
variants were overexpressed in Escherichia coli and 
purified to high homogeneity. The 6xHIS tag was engin- 
eered at the N-terminal of each VirB variants to facilitate 
purification. 

We performed the fluorescence anisotropy experiments 
to study the binding capability of VirB variants to a 
fluorescently labeled DNA FAM-icsB-wt-50 bp derived 
from the ds-acting site upstream icsB promoter (17). We 
observed the binding isotherms that could be fitted to the 
Hill equation (29). The data points of the four VirB 
variants were fitted to the curves with the coefficients of 
determination (R 2 ) ranging from 0.94 to 0.99, respectively, 
indicating the good fittings to the Hill equation 
(Figure IB). VirB FL binds the DNA with an apparent 
dissociation constant K d of 11 ±0.3nM. Lacking the 
N-terminal portion did not affected the DNA-binding 
capabilities significantly, whereas lacking the C-terminal 
portion accounted for > 10-fold loss of DNA-binding 
affinities because VirB delta N binds DNA with a K d of 
12 ± 0.3 nM, and VirB delta C binds DNA with a K d of 
139 ± 15 nM. VirB core harboring a putative HTH motif 
binds the DNA with slightly lower affinity 
(K d =16±3nM) comparing the VirB FL and VirB 
delta N. The aforementioned results demonstrate that 
although the C-terminal portion is important to the 
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Table 1. Data collection and refinement statistics 



Parameter 



VirB core 
(SeMet derivative)/ 
icsb (DNA with 
5-BRU modifications) 
(PDB ID: 3VWB) 



VirB core/ 

icsb (native DNA) 

(PDB ID: 3W3C) 



VirB core/ 

icsp (native DNA) 

(PDB ID: 3W2A) 



Data collection 
Space group 
Cell dimensions 
a, b, c (A) 
oc,P,y (°) 
X ray source o 
Wavelength (A) 
Data range (A) 
Reflections unique 
R sym a (last shell) 
I/oI 

Completeness (%) (last shell) 
Redundancy (last shell) 

Refinement 

Resolution range (A) 
Reflections (non-anomalous), cut-off, 
cross validation 

^work^^free" (last Shell) 

Atoms 

Non-hydrogen protein atoms 

Protein 

DNA 

Solvent 

5-factors average (A 2) 
Protein (A 2 ) 
DNA (A 2 ) 
Solvent (A 2 ) 
r.m.s.d 

Bond lengths (A) 
Bond angles (°) 

Validation 

MolProbity score 

Clashscore, all atoms 

% Poor rotamers 

% residues in favored regions, 

allowed regions, outliers in 

Ramachandran plot 



P2i2!2 

57.32, 163.05, 39.57 

90.00, 90.00, 90.00 

SLS BEAMLINE X06DA 

0.9792 (Se peak) 

46.89-2.42 

27467 

0.067 (0.768) 
21.56 (2.43) 
99.8 (98.8) 
7.77 (7.54) 



46.89-2.42 

27 362 (14781), F>1.99, 1358 
0.2370/0.2710 (0.3640/0.3901) 



2053 
903 
1088 
62 

65.76 
53.25 
77.39 
43.32 

0.004 
0.979 



2.11, 90th percentile" 1 

14.36 

1.90 

96.5, 2.6, 0.9 



¥2^2 

56.14, 163.40, 39.67 

90.00, 90.00, 90.00 

SSRF BEAMLINE BL17U 

1.54001 

40.85-2.43 

13931 

0.049 (0.254) 
20.49 (4.77) 
96.6 (95.8) 
4.95 (3.97) 



40.85-2.43 

13 893 (13 893), F>1.88, 694 
0.2640/0.2797 (0.3532/0.3778 



2061 

931 

1060 

70 

60.1 

48.9 

71.3 

41.0 

0.003 
0.749 



1.76, 98th percentile d 

7.79 

1.90 

97.4, 0.9, 1.75 



C2 

176.32, 39.83, 73.49 

90.00, 90.00, 90.00 

SLS BEAMLINE X06DA 

0.9796 (Se peak) 

48.25-2.78 

23902 

0.090 (0.761) 
12.86 (2.31) 
98.1 (89.2) 
6.68 (5.76) 



48.25-2.78 

23 902 (12 578), F>2.0, 1194 
0.2415/0.2624 (0.4199/0.4428) 



1904 
920 
943 
41 

96.03 
74.13 
118.99 
59.59 

0.004 
0.924 



2.44, 93rd percentile" 1 

18.7 

1.94 

92.9, 6.2, 0.9 



^svm 



where Ihki is the average of symmetry-related observations of a unique reflection. 



X!hk£j |Ihkl,j — Ihkll/J]hklj]jlhkl,j, wnere li lk i 
-^work y~!hkl 

||F obs (hkl)|-|^ calc (hkl)||/Ehkil^obs(hkl)|. 
°Rfree = the cross-validation R factor for 5% of reflections against which the model was not refined. 
d 100th percentile is the best among structures of comparable resolution; 0th percentile is the worst. 



DNA-binding capability of VirB, the N-terminal portion 
of VirB negatively affected the overall DNA binding cap- 
ability (specific and non-specific) of the HTH domain. 
A plausible explanation is that the N-terminal domain 
of VirB might block the DNA binding site of the HTH 
domain in the absence of the C-terminal domain. The 
fitting of the isotherms for DNA binding by VirB FL 
gave the hill coefficient n = 2.8, identifying a positively 
cooperative DNA binding (Figure IB). The removal of 
the N-terminal of VirB seemed to further improve the 
cooperativity in DNA binding (the Hill coefficient for 
VirB delta N is n = 4.0), whereas the removal of C- 
terminal abolished entirely the cooperativity. Because 
VirB delta C and VirB core bind the DNA with the 
Hill coefficients of n = 1.0 and n = 0.9, respectively. 



These results suggest that the C-terminal portion of VirB 
is essential for cooperative DNA binding. 

Crystallization and structure determination 

To gain insights into VirB-DNA interaction, we carried 
out structural studies of various nucleoprotein complexes 
comprising VirB variants and a set of DNAs. We used the 
dynamic light scattering method to evaluate the qualities 
of the samples in solution, identified that only VirB core- 
DNA complex was monodispersed. We next prepared the 
nucleoprotein complexes by mixing VirB core with a 
number of DNA (23-33 bp) containing the c/s-acting 
sites upstream virulence gene promoters. The DNA 
duplexes were blunt-ended or with the complementary 
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Figure 1. The mode DNA binding by VirB variants. (A) VirB variants 
used in this study. (B) Fluorescence anisotropy changes (AA) were 
recorded when titrating VirB FL (black circle), VirB delta N (black 
square), VirB delta C (black down-pointing triangle) or VirB core 
(black up-pointing triangle) into a solution containing fluorescently 
labeled FAM-icsB-wt-50 bp. The experimental points were fitted to 
the Hill equation: A A = AA max *[S] n /(K d n +[S] n ), where AA is the fluor- 
escence anisotropy changes, [S] is the concentration of the VirB variants 
and 'n' is the Hill coefficient. Non-linear regression against cooperative 
binding isotherms reveals dissociation constants (K d ) and Hill coeffi- 
cients (n) of 11 ± 0.3 nM and 2.8 for VirB FL, 12 ± 0.3 nM and 4.0 for 
VirB delta N, 139 ± 15nM and 1.0 for VirB delta C and 16 ± 3nM 
and 0.9 for VirB core. 



sticky ends. A 26 bp DNA 'icsB-wt' and a 31 -bp DNA 
'icsP-wt' (Supplementary Table SI) co-crystallized with 
VirB core. We first determined the crystal structure con- 
taining selenomethionine labeled VirB bound by the 
5-BRU incorporated icsB-wt by single wavelength anom- 
alous dispersion. The resulting model was then used as 
searching model to determine the crystal structures of 
native VirB core/icsB-wt and VirB core/icsP-wt complexes 
by molecular replacement. We could build 116 aa of 122 
aa from VirB core, all 52 bp from icsB-wt and 46 bp of 62- 
bp from icsP-wt in the electron density maps. All crystals 
structures were finally refined to the atomic models with 
excellent refinement and stereochemistry qualities 
(Table 1). 

Overall structure of VirB-DNA complex 

The structures of VirB core from different crystal forms 
are identical. The structural comparison of gave the 
r.m.s.d. values of 0.87-0.95 A over 116 Coc of VirB core. 
In the crystal structure of VirB core/icsB-wt, DNA 
duplexes with complementary 5' overhangs are connected 
end-to-end, forming serpentine helix through the lattice. 
The engineered triple A and triple T overhangs are base 
paired at the junctures (Figure 2A). VirB core molecules 
connect the parallel DNA strands and align along the 



direction perpendicular to the DNA strands, providing 
additional contacts to stabilize the lattice. Similar 
packing was observed in VirB-icsP-wt crystals; however, 
there were insufficient electron densities to build the bases 
at the junctures. 

VirB core comprises seven a-helices folding into the 
N- and C-terminal subdomains. The helices oc2-oc3 form 
a typical HTH DNA binding motif (22). Helices oc4-oc7 
forms a distorted 4-helical bundle that is connected to 
the HTH motif by a linker (Figure 2B). In a search for 
structural relatives of VirB core in protein data bank 
(PDB) using DALI search tool (30), the significant hits 
are ParB-like proteins. VirB core is structurally related 
to the HTH domain of PI plasmid ParB (Supplementary 
Figure SI), F plasmid SopB and RP4 plasmid KorB with 
the r.m.s.d. of 2.0, 2.7 and 4.0 for 113, 106 and 94 corres- 
ponding Coc, (Figure 2B and C). 



VirB preferentially binds the box2 in the inverted repeats 

Both DNA used for co-crystallization contained the 
inverted repeats, boxl and box2 (Figure 3A). The box2 
of icsB-wt icsP-wt are identical (5 / -ATGAAAT-3 / ) that 
match 7/8 positions at the consensus established previ- 
ously (17,24,25). The boxl in icsB-wt (5 / -ATGAAAC-3 / ) 
and the boxl icsP-wt (5 / -CTGAAAT-3 / ) match 6/8 pos- 
itions at the consensus sequence. In the crystal structures, 
VirB core only binds box2 containing with the best match 
to the consensus, whereas the boxl were unoccupied. The 
HTH motif of VirB docks in the major groove of the box2 
in icsB-wt and icsP-wt in the similar way (Figure 3A); 
however, differences were also found. In general, the oc3 
of HTH motif presenting the 'recognition helix' is buried 
in the major groove (Figure 3B) with the helical axis 
running parallelly along the plane of the base. VirB core 
provides six contacts to the backbone phosphates of icsB- 
wt, four of which are clustered on one side of box2, two 
contacts on the other side. VirB core provides seven 
contacts to the backbone phosphates of icsP-wt, four of 
which are clustered on one side of box2, three others on 
the other side (Figure 3A). The asymmetric backbone rec- 
ognition may be essential for the correct docking of VirB. 
It was previously shown that K152 and K164 of VirB are 
essential in DNA binding and transcriptional activation 
(22). Our crystal structure reveals the underlying struc- 
tural basis. N^ of K152 provides a tight salt bridge to 
Op2 of G 8b in icsB-wt (2.61 A), and the backbone NH 
of K152 donates a hydrogen bond to Op2 of the 
adjacent T 7b (2.15 A, 138°). N£ of K164 provides 
another strong salt bridge (2.45 A) to the Op2 atom of 
the phosphate backbone of T 2c . Additionally, we 
identified a number of other residues that also contribute 
to the DNA backbone recognition, including Y151, T166, 
Q170 and F190 (Figure 3A and B). Residues A161 and 
K194 contact the backbone of icsP-wt; however, the cor- 
responding interactions are not observed in VirB core- 
icsB-wt complex. 

VirB provides three base contacts to icsB-wt and two 
base contacts to icsP-wt (Figure 3 A and B). In VirB core- 
icsB-wt complex, the side chain of R167 is oriented to 
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Figure 2. Overall structure of VirB core complexed by the inverted repeats. (A) View of crystal lattice of VirB core/icsB-wt complex. VirB core is 
highlighted in blue ribbon models. (B) Ribbon model of VirB core docking into the box2 of the inverted repeats upstream icsB promoter. Seven 
helices forming the N-terminal HTH motif, a linker and the C-terminal distorted helical bundle are labeled, respectively. The double-stranded DNA 
icsB-wt is shown with stick model with the calculated semitransparent surface, in which box2 is highlighted in orange, boxl is highlighted in green 
and the engineered 5' triple dA is highlighted in black. (C) Structure-based multiple sequence alignment of the transcriptional regulator VirB and the 
plasmid partitioning proteins ParB, KorB and SopB. The secondary structure of VirB is shown on the top, whereas the secondary structure of ParB 
is shown at the bottom. Invariant residues among these proteins are highlighted with red background; conserved residues are shown in red. 



donate two hydrogen bonds to G 3c base (R167 Nr|2 - G 3c 
N7: 1.97A, 167°; R167 Nr|l - G 3c 06: 2.77A, 157°). The 
plane of the guanidinium is twisted with respect to the 
plane of guanine base of G 3c (dihedral angle ~24°) so 
that the side chain of R167 donate another hydrogen 
bond to the adjacent T 4b (R167 Nr|2 - T 4b 04: 2.87 A, 
132°) (Figure 4A) o . R162 donates a hydrogen bond to the 
base of A 7b (2.14 A, 125.5°). The plane of the guanidinium 
is nearly perpendicular to the plane of the adenosine base 
(dihedral angle ~99°) (Figure 4B). In VirB core-icsP-wt 
complex, R167 donates two hydrogen bonds to G 3c base 
(R167 Nr|2 - G 3c N7: 2.10 A, 167°; R167 Nr|l - G 3c 06: 
2.03 A, 160°). The twist angle between the planes of the 
guanidinium and the guanine base reduced to ~14°; thus, 
the side chain of R167 is out of the hydrogen bonding 
range to the T 4b (Figure 4C). R162 adopts a different 
side chain conformation to donate double hydrogen 
bonds to A 7b base (R162 Nr|2 - A 7b N7: 2.70 A, 143°; 
R162 Nr|l - A 7b C8: 2.73 A, 158°) (Figure 4D). Residue 
A 163 provides the third base contact to T 5b in icsB-wt. 
The corresponding contact was unavailable in VirB cor- 
icsP-wt complex. 

The conformational changes of DNA on binding by VirB 

In the crystal structure of VirB core-icsB-wt, there are 
sufficient electron densities to model all 52 base pairs of 



icsB-wt, and the DNA is associated with lower tempera- 
ture factor than that of icsP-wt (Table 1). Therefore, we 
used VirB core/icsB-wt structure to study the conform- 
ation of the bound DNA. The double helix of icsB-wt is 
slightly bent on binding by VirB with the overall curvature 
(31) of ^22°. The bending occurs at the box2 because the 
duplex (C_io-C 0 ) containing an unbound boxl is straight 
(overall bend angle ~5.5°; the average bend of the di- 
nucleotide steps of boxl is 1.4°), whereas the duplex 
(A r Ci 2 ) containing the box2 is bent with the angle of 
~24.0°. The averaged bend of the dinucleotide steps (A 4 - 
T 7 ), an A-tract segment 'AAAT', is 3.3°, in contrast to the 
average bend of 1.5° for the rest base pair steps. In the 
crystal structure of PI ParB-ParS complex, the A2 box 
(5 / -TTG AAAT -3 / ) of ParS shares the most sequence 
homology with box2 of icsB-wt. The A2 box contains 
the identical A-tract segment 'AAAT'. Although the A2 
box is docked by the HTH domain of ParB, it is not bent. 
The average bend of the dinucleotide steps of (A 7 -T 4 ) of 
A2 box is 1.5°. These data suggest that the binding of VirB 
induces a bend at the A-tract segment of box2. The similar 
feature is not shared by ParB. Additionally, the box2 of 
icsB-wt exhibits the evident fluctuation of DNA minor 
groove width that is uncorrelated with appreciable 
changes in major groove width (Figure 5). Comparing 
with the standard B-form DNA, the minor grooves is 
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Figure 3. VirB-DNA contacts. (A) Ladder representations of VirB contacts to 26 bp DNA icsB-wt (left) and to icsP-wt (right) (only the nucleotides 
visible in the crystal structures are drawn). In stands b and c of icsB-wt and icsP-wt, the central base pair is numbered as 0; the base pairs one the 
box2 (orange) side are assigned to positive numbers and the base pairs on the boxl (green) side are assigned to negative numbers. Contacts to 
backbones are indicated with black arrows; contacts to bases are indicated with red arrows. (B) Contacts between HTH domain (blue) of VirB and 
box2 (orange) of icsB-wt. Residues involving the contacts are shown in stick model. Hydrogen bonds to the bases (red) are indicated with dashed 
lines. 



wider at A r G 3 steps followed by the significant compres- 
sion at A4-T7 steps. By contrast, the boxl of icsB-wt that is 
not bound by VirB does not exhibit the similar oscillation 
in minor groove width. In the structure of icsP-wt, we 
observed the similar oscillation in DNA minor groove 
widths in box2 where VirB is bound and in boxl where 
VirB is not bound (Figure 5). However, the boxl of icsP- 
wt is only partially visible in the electron density map, and 
this portion is associated with the highest temperature 
factor. Therefore, the geometry measured for boxl could 
be influenced by its flexibility. 

VirB binds DNA specifically and nonspecifically 

The crystallographic data supports that the box2 of the 
inverted repeats is preferentially recognized by VirB; 
therefore, we investigated the interplay between VirB 
and DNA by EMSA. We used a radioactively labeled 
icsB-wt as the probe and synthesized a mutant DNA in 
which the box2 'ATGAAAT' was replaced by 'GCTGGG 
C, denoted 'icsB-mut' (Supplementary Table SI) as the 
competitor. Such mutations have previously been shown 
to lead to defective VirB-mediated promoter activation 
(17). We also included poly dl-dC in EMSA buffer to 
eliminate non-specific interactions. We observed that 
when the protein concentration increased to 2uM, the 
labeled DNA fully shifted with the VirB variants. VirB 
FL produced a shift of high-ordered complex with large 



size, which could barely enter the native PAGE 
(Figure 7A). However, this large complex did not seem 
to base on the specific DNA-protein interaction because 
the shifted band was efficiently competed with either 
unlabeled icsB-wt or icsB-mut. After the addition of the 
unlabeled DNA, we observed the enrichment of a predom- 
inate shift with smaller size. This band was only competed 
by icsB-wt, but not by icsB-mut, indicating a specific nu- 
cleoprotein complex. These results support our crystallo- 
graphic findings that the box2 of the inverted repeat is the 
key sequence specifically recognized by VirB. Similarly, 
VirB delta N formed the large complex with icsB-wt, 
which could be competed by either icsB-wt or icsB-mut 
(Figure 7B). We observed a lower band after the addition 
of unlabeled DNA, which was competed by icsB-wt, but 
not by icsB-mut. VirB truncations lacking the C-terminal 
portion, VirB delta C and VirB core, did not produce the 
super-shift as observed for VirB FL and VirB delta N, 
suggesting that the C-terminal portion is responsible for 
the formation of the large nucleoprotein complex 
(Figure 7C and D). We observed a single shift produced 
by VirB delta C, which was competed by unlabeled icsB- 
wt but not by icsB-mut (Figure 7C). VirB core produced 
two shifts, suggesting that each box of the inverted repeats 
in icsB-wt may accommodate a HTH domain, although 
with different affinities. The double bands were efficiently 
competed with icsB-wt, but much less efficiently competed 
with icsB-mut (Figure 7D). 
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Figure 4. Base recognition by R167 and R162. (A) R167 reads bases G 3c and T 4b of box2 in icsB-wt. (B) R162 reads bases A 7b of box2 in icsB-wt. 
(C) R167 reads bases G 3c of box2 in icsP-wt. (D) R162 reads bases A 7b of box2 in icsP-wt. Residues (blue) involving the contacts are shown in stick 
model. Hydrogen bonds to the DNA bases (orange) are indicated with dashed lines. 



The determinant for dimerization/oligomerization of VirB 

It was previously proposed that the leucine zipper of VirB 
mediates the dimerization and the C-terminal domain of 
VirB mediates the high-ordered oligomerization (22). Our 
crystallographic study has ruled out the existence of 
leucine zipper domain and showed that the isolated 
HTH domain of VirB exists as monomers in crystal. To 
assess the molecular determinants for dimerization/oligo- 
merization in solution, we carried out the size-exclusion 
experiments to determine the molecular weight of four 
VirB variants. As shown in Supplementary Figure S2, 
VirB variants with the C-terminal domain form dimer/ 
oligomer in solution. The apparent molecular weight 
calculated for VirB FL was 99.2 kDa, which is close to 
the expected molecular weight of a VirB trimer 
(106.2 kDa). The molecular weight calculated for VirB 
delta N was 42.0 kDa, matching the expected molecular 
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Figure 5. Structural distortion of VirB-bound DNA. Major and minor 
groove widths for the icsB-wt (red) and icsP-wt (blue) structures. The 
values are the direct distance between phosphates subtracted by 5.8 A 
for the van der Waals radii of the phosphate groups. Dashed lines 
indicate the groove widths of the standard B-form DNA. 
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Figure 6. Mutagenesis studies. (A) Ribbon diagram of VirB core bound by DNA. Mutations introduced to residues mediating contacts to DNA are 
colored in orange; to control residues not mediating contacts to DNA are colored in green. Leucine repeats with seven-residue intervals are 
highlighted in red. (B) Activation abilities of VirB mutants and truncations. VirB-dependent promoter activation is assessed by measuring the 
P-glactosidase activity of icsP-lacZ fusion. 



weight of a VirB delta N dimer (41.6 kDa). VirB variants 
without the C-terminal domain both exist as monomers in 
solution. The molecular weight calculated for VirB delta C 
was 26.5 kDa, close to the expected molecular weight of a 
monomeric VirB delta C (28.5 kDa); the molecular weight 
calculated for VirB core was 14.0 kDa, matching the 
expected molecular weight of monomeric VirB core 
(14.3 kDa). Collectively, our data support that the C- 
terminal domain of VirB is the determinant for dimer/ 
oligomerization, whereas the N-terminal and the central 
HTH domain do not mediate dimerization. 

Mutagenesis study 

Based on structure, we introduced a set of mutations to 
VirB including mutations at the HTH motif that directly 
contacts DNA and control mutations at the C-terminal 
distorted helical bundle that does not bind DNA. We con- 
structed plasmids expressing VirB truncations, VirB delta 
N, VirB delta C and VirB core. We constructed a VirB 
mutant L203P that is defective in the transcriptional acti- 
vation (22), serving the negative control. We then 
determined the ability of VirB variants and mutants to 



activate the icsP promoter. Plasmids expressing VirB 
mutants were transformed to virB~ S. flexneri 2 a strain 
301 containing an icsP promoter fragment (1256 bp 
upstream of the icsP transcription start site) fused to a 
lacZ gene. We conducted the (3-galactosidase assays as 
illustrated in Figure 6 A and B. VirB truncations lacking 
either N- or C-portion or both failed to activate icsP 
promoter controlling the lacZ gene because the 
measured (3-galactosidase levels fall to background level 
as measured for the negative control L203P. 
Substitutions of the residues involving DNA backbone 
or base recognition, such as mutants Y151A, K152A, 
K152E, A163F, K164A, K164E, T166Y, T166E, Q170A 
and Q170E, lead to the complete losses of VirB-dependent 
activation; however, exceptions were also found. 
Although mutation R162E abolished the (3-galactosidase 
activity completely, mutation R162A had no effect at all, 
indicating that the DNA base contact by R 162 side chain 
is not essential for the activity of VirB. The only essential 
base contact is provided by R 167 because both mutations 
R167A and R167E abolished the activity to the 
background level. Mutation F190E abrogated the 
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VirB-mediated icsP activation, whereas mutation F190A 
had neglectable effect on the activity. This is likely due to 
that the backbone NH group of F 190 mediates the recog- 
nition to DNA backbone by a hydrogen bond that could 
not be disrupted by mutation F190A. Mutation F190E 
may lead to the electrostatic repulsion between the 
glutamic acid and DNA, thus interfering the docking of 
the HTH domain. Control mutations Q176K, K243A at 
the C-terminal distorted helical bundle had little effects on 
the (3-galactosidase activity. 



DISCUSSION 

We determined the crystal structures of the HTH domain 
of VirB complexed by its cognate sites, demonstrating that 
VirB is structurally related to the plasmid partitioning 
proteins. There is no evidence showing that VirB is 
involved in the existing plasmid partition systems (4,32). 
Thus, VirB presents a unique example that a protein with 
ParB-like fold acts exclusively in transcriptional activa- 
tion. It is reasonable to speculate that VirB might be re- 
assigned to a new regulatory role during evolution; 
however, the key structural features of its ancestor are 
preserved. Comparing with the only three available struc- 
tural studies of ParB-ParS-like complexes, the recognition 
of the box2 of icsB-wt by VirB HTH domain resembles the 
recognition of A2-A3 boxes of ParS by ParB HTH 
domain the most. The recognition of the guanosine 
bases (third base from 5' of A-box) by R184 of ParB is 
preserved in VirB-DNA interactions that the correspond- 
ing residue R167 of VirB reads the invariant G 3c . R 162 of 
VirB provides another contact to A 7b base located at 5' of 
box2. However, this contact is not essential for VirB- 
mediated promoter activation. Thus, R167 is the only 
residue conferring the DNA sequence specificity. 

It was previously predicted that VirB processes a 
putative leucine zipper domain (190-230 aa) based on 
the identification of the repeating leucines with 7-residue 
intervals. The putative leucine zipper mediates the dimer/ 
oligmerization, and the substitution of these leucines nega- 
tively affects the dimer/oligomerzation of VirB and its 
activity in vivo (11). Our structures reveal that VirB 
175-246 aa region folds into a distorted helical bundle, a 
structural feature shared with the ParB-like proteins. The 
structure-based sequence alignment indicates that the C- 
terminal distorted helical bundle is less conserved 
compared with that of the HTH motif (Figure 2C), sug- 
gesting the subdomain confers the diversity of function 
among the ParB-like proteins. The structure of VirB 
shows that the repeating leucines are located on the hydro- 
phobic sides of oc5 -oc6 helices with the side chains wrapped 
inside the core of the helical bundle, mediating the hydro- 
phobic contacts among the helices (Figure 6A). Therefore, 
substitutions of these leucines may disrupt the folding of 
the domain. Our current structural data do not support 
that the distorted helical bundle of VirB mediates dimer/ 
oligomerization. Interestingly, the structurally equivalent 
helical bundle in SopB does function as a secondary di- 
merization domain (19), in which oc6-oc7 helices mediate 
the dimer contacts. The oc6 helix in VirB contains a kink 



caused by a proline (P215) in the middle of the helix 
(Figure 2C); thus, the bent oc6 helix cannot support the 
same type of dimerization as seen in SopB, which requires 
a straight helix. Beloin et al. (22) discovered that mutation 
P215A did not affect the activity of VirB and the mutant 
could dimerize more easily than WT VirB in vivo. The 
explanation is likely that the mutation P215A removed 
the bend and straightened the oc6 helix, therefore enabled 
the oc6 -oc7 helices mediated dimerization. 

In the structures of ParB-ParS-like complexes, both 
half-sites of the inverted repeats are occupied by the 
HTH domain (19-21). Intriguingly, the structure of 
VirB-DNA complex show that the HTH domain binds 
only one half site containing the most matches to the con- 
sensus, whereas the other half site is unoccupied. In line 
with this finding, previous studies have demonstrated that 
the complete ParS-like site is not essential for VirB- 
mediated promoter activation, but a close match to the 
consensus is critical (17,23). On binding by VirB, the 
box2 of the inverted repeats exhibits the evident conform- 
ational deviations. The ability of inducing DNA distortion 
presents a distinguished role of VirB, which is not shared 
by other ParB-like proteins because the DNA bound by 
ParB, KorB or SopB usually adopts the standard B-form 
DNA conformation. VirB induces a bend at an A-tract 
segment following the invariant guanine; however, VirB 
does not provide a specific contact to this segment, and the 
base contact provided by R162 has been proven not 
essential for the activity of VirB. It is believed that the 
A-tract segment is highly flexible and may contain a 
static bend in solution (33); thus, our data suggest that 
either the shape or the intrinsic flexibility of the A-tract 
segment rather than its sequence is recognized by VirB. 

The C-terminal domain of VirB shares the least 
homology with the ParB like proteins. Assuming that 
VirB was evolved from a ParB analog for plasmid parti- 
tioning, while the N-terminal and the HTH domains were 
well preserved during the evolution, the C-terminal 
domain might have been diverted to fit its new roles. 
Previous study showed that the C-terminal of VirB 
mediates the oligomerization rather than dimerization in 
a DNA-independent manner, which is essential for its 
function (22). Our bio-SAXS analyses confirmed that 
VirB undergoes concentration dependent oligomerization 
in solution (Supplementary Table S2), which is in accord- 
ance with previously reported cross-linking experiments 
(22). The investigation of oligomerization state of 
various VirB truncations by size-exclusion chromatog- 
raphy also showed that the C-terminal domain is the de- 
terminant for oligomerization. The aforementioned 
experiments were carried out in the absence of DNA, sug- 
gesting that the oligomerization of VirB is independent of 
DNA binding. Thus, the VirB oligomer may be the func- 
tional assembly for DNA binding. In supporting to this 
hypothesis, our fluorescence study showed that VirB 
exhibited a clear sign of positive cooperativity in DNA 
binding that is dependent on the C-terminal domain. 

The EMS A results showed that the C-terminal domain 
of VirB facilitated the assembly of high-ordered VirB- 
DNA complex, forming a super-shift of huge size. 
However, this shifted band was efficiently competed by 
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Figure 7. EMSA of VirB variants. EMSAs of the radio actively labeled icsB-wt (2nM) bound by VirB variants (2uM): (A) VirB FL; (B) VirB Delta 
N; (C) VirB Delta C; (D) VirB Core. Mixtures of VirB variants, labeled icsB-wt and unlabeled competitor (icsB-wt or icsB-mut as specified) were 
incubated and resolved by native PAGE, visualized by autoradiography. Each gel in the figure has the same configuration: lanel, free DNA; lane2, 
VirB variant as indicated; lane 3-8, addition of the increasing amount of unlabeled competitor DNA in 2x, 4x, 8x, 16x, 32 x and 64 x excess of 
proteins. Poly dl-dC was included in each reaction to eliminate the non-specific protein-DNA interactions. 



icsB-mut, indicating that the formation of the large nu- 
cleoprotein complex was mediated by non-specific DNA 
binding. Surprisingly, a smaller band accumulated when 
the competitor DNA were added. The band could not be 
competed with icsB-mut, suggesting a specific VirB-DNA 
complex. Combining with the results from size-exclusion 
chromatography, we speculate that the specific VirB-DNA 
complex might comprise of VirB trimers. VirB core ex- 
hibited the least specificity in DNA binding comparing 
with the other VirB variants, suggesting that isolated 
HTH domain itself only confers partial DNA specificity; 
additional determinants for the recognition of the ex- 
acting site may reside in N- or C-terminal domains 
(Figure 7). 

We demonstrate that VirB binds DNA cooperatively, 
involving both specific and non-specific DNA binding 
with different affinities. Although the C-terminal portion 
of VirB is responsible for the assembly of the high-order 
oligmer, the N-terminal and the central HTH domains 
contribute to the specific and non-specific DNA binding. 
Combining these observations, the DNA-binding mode of 
VirB seems to resemble the model proposed for the DNA 
binding by the leucine-responsive regulatory protein 
(Lrp), a regulator of several E. coli operons (34). 



Lrp binds cooperatively to specific or non-specific DNA. 
The C-terminal domain of Lrp facilitates the oligomeriza- 
tion and is essential for the cooperativity in DNA binding. 
The initial binding of Lrp dimer to DNA induces a bend, 
thus enhances the binding to the next Lrp dimer. 

In conclusion, our results support and extend a previ- 
ously proposed model (16,23) for the relieving of H-NS- 
mediated repression by VirB (Figure 8). The C-terminal 
domain of VirB facilitates the oligomerization, possibly 
the trimerization. VirB oligomer acts as the functional 
assembly for DNA binding. The HTH domain reads the 
c/s-acting site upstream the promoter, mediating the initial 
binding step. The binding by VirB HTH domain induces a 
bend at the A-tract segment, which assists the subsequent 
DNA-binding steps, initiating the wrapping of DNA 
around VirB oligomer. The intrinsic flexibility of the 
invariant A-tract is likely important to the wrapping of 
DNA. After the initial binding at the ds-acting site, the 
subsequent DNA binding is mediated by non-specific 
VirB-DNA interactions. The docking of the first VirB 
oligomer at the ds-acting site serves as a nucleation site 
that subsequently recruits more VirB oligomers. Thus, 
VirB polymerizes alone DNA strand (27), forming the 
large non-specific VirB-DNA complex. The wrapping of 
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Figure 8. Proposed mechanism for the alleviating of the H-NS 
mediated transcriptional repression by VirB. The C-terminal domain 
(yellow) mediates the formation of VirB oligomer (blue). The trimer 
drawn here represents a given oligomeric form of VirB in solution, 
possibly trimeric. The HTH domain of VirB recognizes a c/s-acting 
site upstream the promoter repressed by H-NS-DNA complex. The 
binding by the HTH domain induces distortion and bends DNA to 
initiate DNA wrapping around VirB oligomer. The docking of the 
first VirB oligomer at the ds-acting site serves as a nucleation site 
that subsequently recruits more VirB oligomer; thus, VirB polymerizes 
alone the DNA. The DNA remodeling by VirB followed by the 
spreading of VirB disrupts DNA-N-HS-DNA bridges, leading to the 
release of the trapped RNA polymerase (orange). The transcription of 
the virulence gene is de-repressed. 



DNA and the polymerization of VirB eventually destabil- 
izes the H-NS-DNA complexes, leading to the collapse of 
the repressive H-NS-DNA complexes so that the trapped 
RNA polymerase is released, and the transcription of the 
virulence genes resumes. 
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